PM/98-08 
THES-TP 98/04 
hep-ph/ 9803422] 
March 1998 



Testing the Higgs boson gluonic couplings at LHCfl 



G.J. Gounaris", J. Layssac 6 and F.M. Renard & 

department of Theoretical Physics, University of Thessaloniki, 
Gr-54006, Thessaloniki, Greece. 

^Physique Mathematique et Theorique, UMR 5825 
Universite Montpellier II, F-34095 Montpellier Cedex 5. 



Abstract 

We study Higgs + jet production at hadron colliders in order to look for new physics 
residual effects possibly described by the dim = 6 operators Ogg an d Oqg which in- 
duce anomalous Hgg and Hggg couplings. Two ways for constraining these operators 
at LHC may be useful. The first is based on the total Higgs boson production rate 
induced by gluon-gluon fusion, in which the main cause of limitations are due to theo- 
retical uncertainties leading to sensitivities of \dc\ — 3. x 10~ 4 and \da\ — 1.4 x 1CT 3 for 
the corresponding anomalous couplings, in the mass range 100 GeV < m# < 2 OOGeV. 
These results imply sensitivity to new physics scales of 51 and 24 TeV respectively. The 
second way investigated here concerns the shape of the Higgs transverse momentum; for 
which the theoretical uncertainties are less severe and the limitations are mainly induced 
by statistics. A simple analysis, based on the ratio of the number of events at large and 
low pt at LHC, leads to similar sensitivities, if only the H — > 77 decay mode is used. But 
the sensitivities can now be improved by a factor 2 to 10, depending on the Higgs mass, 
if the Higgs decay modes to WW*, ZZ*, WW, ZZ are also used. 
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1 Introduction 



The Higgs mechanism is the cornerstone of the Standard Model (SM) [|IJ. Much experi- 
mental and theoretical effort has already been done to discover the Higgs boson, and to 
understand the origin of the scalar sector and constrain the Higgs boson mass, which is 
to a large extent arbitrary. At present, the negative experimental searches at LEP2 lead 
to the conclusion that m# > 77 GeV 0, while precision measurements at Z peak give 
mu = 115^gg 6 Ge\^ within the SM context These values are well within the bounds 
obtained from purely theoretical considerations ||. 

Searches are presently going on at LEP2 and the Tevatron and will be pursued at LHC. 
If the Higgs boson is light, its discovery should be at hand. However, such a discovery, 
will only constitute the first step in the study of the scalar sector. Detailed studies of the 
properties of the Higgs candidate should then be performed aiming at first confirming the 
Higgs nature, and then trying to get information about the origin of the scalar sector. 
This origin may very well lie in an underlying fundamental dynamics generically called 
new physics (NP). 

Assuming that the NP degrees of freedom are associated to an energy scale Ajvp much 
larger than any conceivably accessible energy range, any observable NP manifestation 
could only be in the form of residual effects modifying the SM interactions among usual 
particles and the (yet to be discovered) Higgs boson. At the energy range of the foreseeable 
colliders, the leading such effects may arise from the set of dim = 6 SU(3) x SU(2) x U(l) 
gauge invariant operators [El, |5], [7|, involving an isodoublet Higgs field and the particles 
with the highest affinity to it, which inspired by SM, are taken to be the quarks of the 
third family. Of course, apart from these fields, the NP operators necessarily also involve 
gauge bosons, inevitably introduced by the gauge principle whenever a derivative appears. 
The complete list of the CP conserving such operators has been given in 0, , while the 
CP violating ones have appeared in || |[|. Tests of their effects at present and future 
colliders have been given for most of them in || |7|, [K| 11, |l^, |13|] . 

Among these operators, those inducing anomalous Higgs-gluon couplings beyond the 
SM ones arise when the NP degrees of freedom are coloured. Examples of dynamical 
models in which such operators are generated have been pointed out ||, but very little 
has yet been said about their possible signatures fll4f . 

The aim of this paper is to study the Higgs-gluon couplings induced by these operators, 
namely the CP-conserving one dubbed Ogg and its CP- violating analogue Ogg, respec- 
tively associated to dimensionless couplings do andjic- These operators are the gluonic 
analogues of the operators Obb, Oww, Obb and Oww, which are inducing anomalous 



Higgs couplings to the electroweak gauge bosons and have been studied in |T^, [16], [L7|, §]. 
Since the SM couplings among Higgs bosons and gluons, first appear only at the 1-loop 
level, (essentially through a top quark loop), and are therefore somewhat reduced, there 
is a favorable situation for the detection of any NP contribution induced at tree level by 
the possible operators Oqg and Ogg- Other types of NP effects, like e.g. the anomalous 
gtt and Htt couplings modifying the H gg and Hggg interactions through a top-loop con- 
tribution, may also appear, 113]. Below, we show how the presence of these NP effects 
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can be tested through studies of Higgs boson production and decay. 

As the H — > gg branching ratio only represents a small fraction (6-7% for inn ~ 
100 — 150 GeV and much less for higher masses), accurate measurements of the Hgg 
couplings through decay width T(H — > gg) will need a copious Higgs boson production 
and very powerful methods for distinguishing gg from light quark final states. At e + e~ 
colliders, only a few H —* gg events are expected through HZ production; i.e. a few tens 
at LEP2 and a few hundreds at a higher energy linear collider (LC), depending on the 
achieved luminosity. An analysis based on a luminosity of 50 at a 500 GeV LC, 
led to an uncertainty of ±39% for the sum cc + gg (Tj|. The same situation probably 
also arises for H production in a 77 Collider. At the upgraded Tevatron, a hundred of 
Higgs particles should be produced (mainly through the process pp — > WH + X) for 
mn — 100 — 150; but its observation seems possible only in the bb mode |20| . 

The situation should be quite different at LHC where a copious Higgs production of 
about 10 5 - 10 6 events per year for m H < 200GeV, should be expected |22], |23|, |2§ [25|. 
At the LHC energies, the largest cross sections arise from subprocesses with gg, WW/ZZ 
and qq initial states, as well as from bremsstrahlung off a top quark. Among these, the 
process gg —* H largely dominates. So this should be the best place to look for anomalous 
Hgg, Hggg couplings. The standard prediction for this process has been computed and 
found very sensitive to higher order QCD corrections due to soft gluon effects, for the 
resummation of which there exist large theoretical uncertainties [p6j |27| . This means 
that an accurate measurement of Hgg couplings through the total production rate is 
only hindered by theoretical uncertainties, while statistics are enormous. Assuming a 
conservative theoretical error for the QCD corrections, we estimate the discovery limits 
for the Ogg and Ogg effects. These limits are interesting, since they correspond to NP 
scales lying in the tens of TeV range. 

A theoretically cleaner (and in any case complementary) study of the operators Ogg 
and Ogg could be achieved by looking at the Higgs transverse momentum distribution 
in the Higgs + jet production process. To our knowledge, this has never been discussed 
before for the search of NP effects. To achieve this, one first needs the 1-loop SM am- 
plitudes for gg — ► H g, gq — > Hq and qq — ► Hg, that will interfere with the NP ones. 
These SM contributions have been computed in p5|, E9l. We have recomputed them and 



checked numerically the agreement with the results obtained previously. Adding then the 
NP contributions due to the operators Ogg an d Oqg, we have examined how the SM 
predictions for the various observables, (like e.g. the distributions of the Higgs rapidity 
and transverse momentum, and the Higgs+jet invariant mass and angular distributions) 
are influenced by NP. 

The most striking effect is, as expected, in the pr dependence. The SM prediction (as 
well as the one of any other NP model affecting only gtt and Hit couplings) should drop 
down as soon as pt > m t ; whereas the NP contribution due to Ogg an d Ogg operators 
(associated to a heavy A^vp) stays flat, leading to a clear signal for anomalous couplings. 
We then derive the observability limits on do and cIg by considering the ratio of Higgs+jet 
production rates, at the high and low pt regions. 

The content of the paper is the following. In Section 2 we present the Ogg and Ogg 
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operators and derive the unitarity constraints relating the NP couplings to the scale A^p, 
where the new degrees of freedom start being excited. We then present the formulae 
for the NP and SM contributions to the H — > gg decay, and to the Higgs boson and 
H + jet production in pp collisions. The sensitivity limits to the NP couplings are 
derived in Section 3, first on the basis of the expected accuracy in the measurement of 
the T(H — > gg) width; and then on the basis of the Higgs production rate and the various 
H+ jet distributions at LHC. We find that for 100 < m H < 200GeV, the study of the total 
Higgs production rate, and the study of the p^-distribution in the case that only H — > 77 
is used for if-detection, give comparable results. On the other hand the p^-distribution 
technique starts becoming superior whenever the H — > WW*, ZZ* , WW, ZZ decay 
modes are also used to increase its statistics, which may even lead to an order of magnitude 
improvement as ran approaches the 200 GeV region. Section 4 summarizes the results 
and their implications for the NP search. Technical details on invariant amplitudes, loop 
computations and parton kinematics are collected in an appendix. 

2 Formalism. 

2.1 The NP operators O gg and O gg . 

We consider the NP effects arising from the effective Lagrangian 

C-np = —^-Ogg + -^-Ogg > (1) 
where dc and dc are NP dimensionless couplings, and the NP operators are 

Oca = (V*-^)^^ , (2) 

Ogg = (*'*) G„ F ■ '3"" , (3) 
with G%, = l/2e„„ A „G"" and 

* = ■ (4) 

Unitarity allows to establish for each operator an unambiguous relation between the 
NP coupling constants appearing in ([!]), and the corresponding energy scale A^p, at which 
unitarity is saturated. This scale supplies a practical definition of the scale where the new 
physics generates the corresponding operator || [30|, |16[ . 

In the case of the Ogg and Ogg operators, the strongest unitarity constraint arises 
from the J = partial wave transition amplitude affecting the colour singlet channels 
\gg ± ±) and \HH). Using ([I]), we thus get for Ogg the unitarity relation 



with c = 31 for da > 0, and c = — 1 for d G < 0. The unitarity relation for Oqg is 

47T / V 2 X 



2.2 The i7 — > gg decay width 

In this section we give the modification to the Higgs gluonic decay width induced by tree 
level effects of the operators Oqg an d Oqg- 

The Hgg gauge-invariant amplitude^ R(Hgg) is given by 

R(Hgg) = if(k 2 )5 ab [(k.k')(e.e')-(e.k')(e'.k)} , (7) 

where the gluon momenta and polarization vectors are denoted as (k, e) and (k 1 , e') 
respectively^. In (^), the momentum k is allowed to be off-shell, while the other gluon 
with momentum k' and polarization e' is always on-shell. Under such conditions, there is 
only one gauge invariant form for Hgg, which is shown in (|7p. The indices (a, b) specify 
the colours of the two gluons. In SM, the dominant contribution arises from the top 
triangle in Fig. la and is given by 

f{k 2 )=^m) , (8) 

TTV 

in terms of f(k 2 ) presented in Appendix Al. Here we only note that /(0) agrees with the 
result quoted in []23], ^4], [26|, ^8|. In the heavy quark limit m t ^> m#, this leads to the 



well known result [31, 29 



1 



/(°)-"^ ■ ( 9 ) 



3m 



In the presence of the NP contribution given in (JJ), the Higgs decay width into 2 real 
gluons is given by 

r(^ w ) = g{| /( o)-^ + $b|'} . (io, 

As expected, only the CP-conserving NP contribution interferes with the SM one. 

2.3 Higgs and Higgs+jet production at pp colliders 
2.3.1 Single Higgs production 

At lowest QCD order, including the effect of the NP operators Oqg and Ogg, the total 
rate for pp — * H + X due to the subprocess gg — > H is given by 

a (pp^H + X) = a r H L g9 , (11) 



1 The phase of the amplitude is defined to be that of the S-matrix element. 
2 A11 momenta are taken as incoming. 
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where 



J 99 



1 dx 



*(x)»(2) , 

X \ X J 



Ti- 



ll 



a 



6m 



-T(H^gg) = -{\f(0) 



H 



4da 2 16- 
1 +^tI«g| 



(12) 
(13) 



t# = m 2 H /s and g(x) is the gluon distribution function inside the proton. 

QCD corrections correspond to including loop corrections to gg — > H, as well as 
corrections due to associated production of massless partons together with the Higgs, in 
the processes gg — > H g, qq — ■> Hg, gq — ► f/g, gg — > Hq [^T], |22], [!2|, ^4], |2?|, |26], |33|. The 
result in SM is 



a(pp -> # + X) = cr c 



7T . 



t^L™ + AcTgo + Acx w + Acr w - + Act, 



'gq 



(14) 



where the various terms correspond to gg, gq(q) and qq initial state. Depending on the 
value of the Higgs mass, these QCD corrections increase the Higgs production rate by 
60% to 90% Km 02, 123, m E7L R3| . 



2.3.2 H+jet production at large pr- 



We now turn to Higgs + jet production (Fig.2), which at a hadron collider, takes place 
through the subprocesses (Fig. 3) gg — > Hg, gq — > Hq, gq — > Hq, qq — > Hg |28|, j29fl . 

We have repeated the computation of of the triangle and box contributions to the 
various types of subprocesses participating in the H + jet production (compare Fig. 3), 
and we have added the tree level NP contributions due to Oqg and Oqg- The Oqg 
contribution interferes with SM, whereas the Oqg one does not and adds quadratically 
in the cross section. 

Details on the amplitudes are given in Appendix A2; while those on the kinematics 
of the two-body inclusive parton model distributions, are presented in Appendix A3. We 
collect there the expressions for the Higgs transverse momentum distribution da/dpr, the 
H rapidity distribution da/dyn, the (H+jet) invariant mass distribution da/dM, and the 
angular distribution da/d\- They are obtained by convoluting the elementary differential 
cross section da/dt for the various subprocesses, with gluon and quark distributions. 

For gg — > Hg the elementary differential cross section is 



dojgg -» Hg) 
dt 



3ut 



-us 



A^+UAA 2 



222\ 



A(\A 2 \ 2 u 2 + \A s \ 2 t 



AuiU{A 2 Al) + 2uM{{tA 3 - uA 2 ){A\ + sA\)) + siiuf^A^l) 



+ 



3a s |d 



G\ 



2p 



_ In _ Mky 

s s 



+ 



Pt 



(1 



+ 1 + 



,m 



H\4 



(15) 



In this expression the functions Aj contain both SM and Oqg contributions as given by 
( [A. 18 , A.19 , A~T3" , A.14 ) in the Appendix A2. The quantity pt denotes the transverse Higgs 
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momentum; compare Appendix A3. As expected, the CP-violating Oac contributions 
does not interfere with SM and appear separately in quadratic form. This da quadratic 
term is the same as the quadratic term of the Oac contribution (coming from the functions 
Ai) by just replacing d G — > da- 

While for qq — > Hg and gq(q) — > Hq(q) we get 

d&(qq — > Hg) a s ~ 2 -2\fi./v-\ 4g?g,2 16. 



(i 2 +u 2 ){\m-—\ 2 +^\dc?\ , (i6) 

I v v z J 

§ 2 + u 2 )\\m-^\ 2 + 1 4\da\ 2 } - (IT) 



di 36s 3 

da(gq — > Hq) da(gq — > Hq) ot s 2 -^fir/^ 4g?g, 2 16, ~ 2 



dt dt 96s 2 |t 

The corresponding expressions for gg — > if(yf, g(7 — > Hq, qg — > Hq are obtained through the 
replacement t <^=^> m. The function f(k 2 ) representing the Standard Model Hgg triangle 
loop contribution, has been defined in @ and in Appendix Al. 

3 Sensitivity to NP couplings 

3.1 From T(H — > gg) measured at Linear Collider. 

We give an estimate of the sensitivity to the NP couplings da and da for a light Higgs 
boson (ms ~ 100 — 200^6^), assuming a global uncertainty on T(H — > gg) of about 40%. 
This would cover the theoretical uncertainties on a s and higher QCD effects, as well as 
the experimental errors in the measurement of this decay width. It seems that even with 
a high energy, high luminosity e + e~ Linear Collider (LC), this is the best one can expect 



[IE]. With this assumption one obtains the sensitivity limits 

1^1=^1/(0)1 , Mg| 2 = ^|/(0)| 2 (18) 
which, using also the unitarity constraints (|5|, ||), means f\ 
\d G \ = 6. x 10~ 4 (A NP = 36 TeV) , |d G | = 2. x 10~ 3 (A NP = 20 TeV) . (19) 

These would be already quite remarkable values. They arise because one compares 
tree level NP effects with 1-loop SM contributions. 



3.2 From the Higgs production rate at LHC. 

At LHC, Higgs production should be dominated by the gg — > H process. We therefore 
consider the production rate in pp — > H + X given by ( |TT| - |TJ]). Theoretically, the 
sensitivity to NP should be similar to the one expected from T(H — > gg), as it is this 

3 Note that our notation for the NP scale is different from the one in Il3]. 
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same quantity which controls the production rate. Experimentally, this Higgs production 
process is very interesting, due to the large number of events expected in the various Higgs 
decay modes. For a light Higgs (m# ~ 100 — 150 GeV) the 77 channel is experimentally 
favored (since the statistics is enormous anyway), and the expected experimental accuracy 
is of a few percent |20| . 



In this analysis we assume that the NP effect on Br(H — > 77) is negligible. This 
stems from the observation that it could only come from two possible sources. The first 
one is a direct NP effect on the #77 coupling arising from other NP operators (not 
involving gluon fields), which can be constrained through different processes though, as 
discussed in |15|, [L6|, [L7|, || . The second possible source comes from a modification of the 
T(H — > gg) width induced by O gg , @gg, which in turn changes accordingly the total 
if- width and Br(H — > 77). However, as the Br(H — > gg) is only a few percent, even 
a 20% modification of T(H — > gg), induces only a ~ 1% variation to Br(H — > 77), 
which is negligible compared to the experimental uncertainty. In addition, to the already 
mentioned theoretical uncertainty in the QCD corrections, one should add the uncertainty 
in the parton (mainly the gluon) distribution functions. In [ 26 , B3 , [J5| , these are estimated 



to be of the order of ±10 % to ±20 %. So they turn out to be the dominant source 
of uncertainty in this analysis. Assuming a global uncertainty of ±20 % will give the 
sensitivities: 

\d G \ ~ 3. x 10~ 4 (51 TeV) , \d G \ ~ 1.4 x 10~ 3 (24 TeV) (20) 

Numerically, these numbers are only slightly better than the ones given in the preceding 
Section on the basis of T(H — > gg) at LC. But since LHC should anyway run before LC, 
we must emphasize the importance of these results, which give indeed very interesting 
constraints on the NP effects. 

In concluding this subsection we note that the main difficulty here comes from the- 
oretical uncertainties, which prohibit us to make a better benefit of the huge statistics 
brought in by LHC. This is why in the next Section, we turn to a theoretically cleaner 
NP signal, free of the normalization uncertainties of the total cross section. 



3.3 From the shape of the H+jet distributions at LHC. 

We now look at the NP signal based on the relative (but drastic) differences in the shape 
of the H production at large transverse momentum. 

At Pt ^ m t, the SM distribution (based on the triangle H gg and the box Hggg 
contributions due to the top quark loop in Fig.l) starts being sensitive to the non- 
locality of the Hgg vertex and falls off; whereas the NP contribution (associated to a 
large scale A^p) is still local and remains flat. In fact, one can check that the pure NP 
contribution behaves in the same way as the contribution called SM e ff in Figs. 4-7 below, 
which is derived by taking the large m t limit in SM and identifying \d G \ = or \d G \ = jrf-, 

The results for various (H + jet) distributions are shown in Figs. 4-7. In all our 
illustrations m# = 100 GeV is used, but we have checked that there is very little change 
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when this mass increases up to 200 GeV. In Figs. 4-7, we compare the complete 1-loop SM 
prediction (labeled SM), with the large m t approximation to SM (labeled SM e ff), and 
the effects due to the presence of the O gg operator with d G = +10~ 3 (labeled +d G ) or 
d G = — 10~ 3 (labeled —d G ), as well as the effect of the O gg contribution for \d G \ = 10~ 3 
(labeled d G ). The results include all possible subprocesses due to gg, gq gq and qq initial 
states; compare ( |A.45| - [A"^48l) . As shown in Fig.6, the subprocess gg — > Hg dominates the 



Pt distribution for x T = 2E TH /^J~s < 0.1; while for x T > 0.1 it is gq,gq that dominate; 
the qq contribution being always smaller. 

The main features concerning the NP observability are discussed below: 

Shape of yu and x distributions. 

The shape of the yn and x (compare ( |A.42| ) distributions in Figs. 4, 5) do not seem to be 



notably different in the SM and in the NP cases. They only differ in absolute magnitude, 
roughly in the same way as the total production rate. So, from these measurements, we 
cannot expect an improvement in the determination of the NP couplings as compared to 
the one obtained from the total production rate. 

Shape of xt and M distributions. 

On the contrary, the transverse energy and the H + jet invariant mass distributions 
(Fig.6, 7) are very sensitive to NP. At large xt or M, the NP contributions differ from 
SM by a flattening of the distributions. As already stated, such behaviour is due to the 
locality of the NP interaction implied by its high A^p scale. For the same reason, this 
is also true for the SM e ff approximation describing the large m t limit of the SM, which 
should thus become inadequate for p T > m t . The size of the NP effect, is fixed by the 
value of the d G and d G couplings, which once they are determined also fix the accessible 
values of Atvp through the unitarity relations (|, |). 

We also remark that the quadratic O gg contribution always increases the rate, whereas 
the linear O gg one produces a constructive or destructive interference with SM, depend- 
ing on the sign of d G . The observation of a destructive effect would be a clear indication 
for O gg . 

Ratios and sensitivity to NP 

As already stated, this sensitivity is based on the change in the shape of the pt distribution 
induced by NP. In order to quantify it, we consider the ratio 

R = N High /N Low , (21) 

where Nuigh and Niow are the number of (H + jet) events in the high and the low trans- 
verse energy domains [xi, XTmax] and [xTmin, , respectively; compare Fig.6. The lowest 
and highest points XTmin, XTmax are chosen as xxmin = 0.0257 and xxmax = 0.25; but the 
results are independent of their precise values. Note in particular that x-rmm = 0.0257 
corresponds to pt — 150 GeV, for which the use of the leading QCD expression should 
be adequate. The intermediate value is chosen in order to maximize the sensitivity of 



9 



R to NP effects. It naturally turns out that Xi lies close to the value where the SM e ff 
prediction crosses the exact SM one. We choose therefore Xi = 0.05. 

As seen in Fig. 6, the number of events in the low xt domain is determined solely 
by the SM contribution; while the NP effects mainly influence the events in the high xt 
domain. This is even more pronounced for xt > 0.25, but there are very few events there 
anyway. The NP observability limit is then defined by demanding that the NP effect on 
the ratio R is larger than the la statistical fluctuations of it. 

At LHC, with an integrated luminosity L = 100 fb' 1 per year, and 3 years of running 
for the two experiments ATLAS and CMS, we find N%% h = 22320 , Nf^ = 417000 at 
the two Xt domains. This should a priori lead to an excellent sensitivity to NP. However 
these numbers are reduced by the branching ratios of the Higgs boson to the observable 
channels and by detection efficiencies. For the mass range 100 GeV < mn < 200 GeV, 
the relevant subprocesses are gg — > H — > 77, gg — > H — > WW* and gg — > H — > ZZ* [20 . 
If we only use the H —>■ 77 mode, whose branching ratio is very low (about 2 x 10~ 3 for 
m H ;$ 150 GeV and negligible for higher masses), we obtain Nfjf gh = 22 , = 420. 

With such numbers the statistical uncertainty on the ratio R is about 25%, which leads 



to sensitivity limits for the NP couplings similar to those given in (|20|). 

This result can be substantially improved for most of the mass range 100 GeV < 
m H % 200 GeV, by using also the Higgs decay modes to WW*, ZZ* which strongly 
increase with m#. Thus the WW* {ZZ*) modes reach e.g. at m# ~ 150 GeV the level of 
50% (5%) respectively. Above the corresponding thresholds, the WW and ZZ modes will 



of course dominate and greatly increase the statistics. In [^0| an analysis of the two-weak 



boson final states is performed, which shows that their statistical error is already smaller 
than the one for the 77 case, in most of the above mass range. Considerable work should 
of course still be done on the detection of the various Higgs decay modes, before we fully 
identify all potentialities 0, Using though the existing results we are led to the 
conclusion that an improvement of the sensitivities given in (^) by a factor of 2 to 10, 
depending on should be possible. This is particularly true for m# > 120GeV. 

We should also note that our present study of the pt shape through the ratio R, is 
intended only as a preliminary orientation. In the actual search, a global study of the 
shape of the pr spectrum should be done in a more precise way during the event selection, 
taking into account the pt dependence of the background and all characteristics of the 
detectors Q. This is beyond the scope of this paper and our competence. We nevertheless 
believe that our simple study has shown that it is reasonable to expect an appreciable 
improvement in the sensitivity to the NP contribution. 

Finally we note that we obtain similar results by applying the same procedure to the 
invariant mass distribution. The intermediate value which optimizes the sensitivity is now 
found as Mj ~ 0.9 TeV. The results for do and do sensitivities turn out be similar to the 
ones obtained from the xt distribution. 



4 A recent study|37| of the process pp — * 77 + jet shows that the background has indeed different 
distributions. 
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4 Conclusion 



In this paper we have studied the dim = 6 SU(3) x 577(2) x U(l) gauge invariant operators 
O gg an d O gg , which induce anomalous Higgs boson gluon interactions. Such effects 
should appear when the NP degrees of freedom carry colour and simultaneously couple to 
the Higgs sector. For studying the anomalous Hgg and Hggg couplings, we have looked 
at the tests which can be realized from Higgs production and decay. 

At e + e~ colliders, the measurement of the Hgg coupling has to rely on the difficult 
measurement of the 77 — > gg width. At LEP2 only a few events could be observed. At a 
500 GeV LC, for uih — 100 — 150 GeV, with a few hundreds of events accessible, one can 
only expect to reach the sensitivity limit 

\d G \ ~ 6. x 10~ 4 (A NP ~ 36 TeV) , \d G \ ~ 2. x 10~ 3 (A NP ~ 20 TeV) . (22) 

At LHC, the ggH coupling controls the 77 production rate. Statistics is huge. The 
main problem in this case is due to the large uncertainties (about 20%) in the QCD 
corrections to the total cross section. This limits the sensitivity to NP to the values: 

\d G \ ~ 3. x 10" 4 (A NP ~ 36 TeV) , \d G \ ~ 1.4 x 10~ 3 (A NP ~ 20 TeV) , (23) 

for 100 < ran J$ 2QQGeV . The importance of this LHC measurement, which in fact 
should precede the aforementioned one at a Linear Collider, should be appreciated. 

The sensitivity to the NP couplings will be notably improved if at LHC we also study 
the shape of the 77 + jet distributions at large transverse momentum or large invariant 
masses. The change which is induced in the shape of these distributions constitutes a clear 
signal for NP. We have explored this possibility by considering the ratio of the number 
of events at large transverse momentum or invariant mass to the one at low values. This 
quantity is rather insensitive to theoretical uncertainties (due to SM itself or to anomalous 
gtt or Htt couplings), but the precision is now controlled by the smallness of the expected 
statistics. If only the 77 — ► 77 mode is used, we reach about the same sensitivity as the 
one given in ( p3|) . This is already a very interesting result as it constitutes a theoretically 
clean independent measurement. But now, depending of the value of the Higgs mass, 
the sensitivity can be notably improved by using the WW* and ZZ* modes, as well as 
the real WW and ZZ modes as soon as m# goes above the corresponding thresholds. 
Depending on the value of the Higgs mass, we then expect an improvement by a factor 2 
to 10 in these observability limits. 

The values of the NP scales that such an analysis should allow to reach i. e. 50 to 100 
TeV, are remarkable. They are of the same order of magnitude as the ones expected in the 
electroweak sector for the operators Oww, @ww, Obb and Obb describing anomalous 
HZZ, HZ'-f and ^77 couplings and which affect the corresponding H decay modes as well 
as HZ, #7 production in e + e~ collisions and 77 production in 77 collisions, [|, [1^, [16], [L7|| . 

A comparison of the effects or of the limits obtained in these two (electroweak and 
gluonic) sectors should tell us about the flavour and colour content of new physics. 
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Appendix Al: 1-loop SM contribution to the ggH amplitude 



The form factor f(k 2 ) determining, (through the diagrams in Fig. la), the Hgg coupling 
in SM in the case that the momentum k of one of the gluons may be off-shell, is 



f(k 2 ) = C {k 2 ) 



Am 2 



m 2 H — k 2 



+ ■ 



2k 2 



,F(k 2 



(A.I) 



m 2 H — k 2 {iTi 2 H — k 2 ) 2 
where all conventions are given in Sect. 2. 2. The triangle loop is computed through the 



FF-a package |38| using the standard Passarino-Veltman method |39| and the notations 

C (k 2 ) = C (k 2 , 0, m 2 H ; m t , m t , m t ) . 

[(1 — x)m\ + xm\ — x{l — x)m 2 H — ie] 



of @ and gTj with 
and 



(A.2) 



F{k z ) = In 



[(1 — x)m\ + xm 2 — x(l — x)k 2 — ie] 



(A.3) 



In the limit k 2 — > 0, in which both gluons are on-shell, f(0) matches the result 



7(0) 



2m 2 



(1 + (1 - r t )f(r t )) , 



(A.4) 



where r t = 4m 2 /m 2 H and 



f{T t )=[sin- 1 {l/^n)\ 



Kn) 



In 



1 + 



n 



l-VT 



m 



if 



if 



T t >l 



T t <l 



(A.5) 



Appendix A2 : CP conserving amplitudes for gg — ► Hg 

Below we discuss the 1-loop SM contribution, as well as the CP conserving NP con- 
tribution, to the invariant amplitudes for the process 

g(e, k, a) + g(e, k', b) -> H(q) + g(e, q', c) , (A.6) 

where the polarization, momenta and colour indices of the gluons are indicated in paren- 
theses. The 1-loop SM contribution to this amplitude has been computed long ago by 
R.K. Ellis et.al. |p9| . It arises from the triangle and box terms appearing in Fig.l and 
involved in Fig. 3. We have recomputed their result in order to make sure that the NP 
contribution is added correctly to the SM one. 

If Bose-symmetry among the three external gluons were ignored, there would had been 
14 different Lorentz invariant and CP conserving forms contributing to the amplitude, 
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which are reduced to just 4 ones, when gauge symmetry is imposed. For the gluon 
momenta and polarizations indicated in ( |A.6| ), these are 



N, = [(e.k)(k f .q') - {e.k'){k.q')\{{k.q%e.e'){k' .q') ~ (e.k')(e' .</)] 
He'.q')(e.q')(k.k')-(e.q')(e'.k)(q'.k')} , 

N 2 = [(e.e)(k.q') - (e.q')(e.k)W .k)(k' .q>) - (e' .q')(k.k')} , 
N 3 = [(e'.e)(k'.q f ) - (e' .q')(e.k')][(e.k')(k.q') ~ (e.q')(k.k')] , 

N 4 = [(e.k)(k'.q') - (e.k')(q'.k)][(e.k')(k.q') - {e.q'){k.k')\ ■ 
[(e'.k)(k'.q')-(e'.q f )(k.k')\ . 

As a result the total contribution to the amplitude is written asQ 

4 

R(gg -> Hg) = -fabc ^2 a j(^ u) N i > 

where 



(k + k'f, t=(q-k) 2 , u = (q'-kf 



(A.7) 

(A.8) 
(A.9) 

(A.10) 
(A.ll) 
(A.12) 



with the momenta defined in ( |A.6| ). The Aj amplitudes in ( [A.ll] ) are not the same as the 
ones used in |p9|| . In accordance with this reference, the requirements of Bose-Einstein 
statistics among the three external gluons reduce the number of these amplitudes to two, 
by determining A%, A3 in terms of A\, through 



A 2 (s,t,u) 
A 3 (s,i,u) 



ts 



A x (u,t, s) 



and imposing the constraints 



Ai(s,t,u) 
A 4 (s, i, u) 

A 4 (s, i, u) — A 4 (t, s, u) 



us „ A . 



Ai(s,w,f) , 
A 4 (s,u,i) , 

2 2 

-Axft, s,w) Ai(s,t,u) 

t s 



(A.13) 
(A.14) 



(A.15) 
(A.16) 

(A.17) 



Thus, we only need to give the 1-loop SM and the tree level Oqg contributions to A\ 
and A4, which are 



A x (s,i, u) 



7t 2 vtu 2 



t — rriir-,,., u-,^, 1- 



32g s d G \ ( u + s)(i+s) 
st 2 u 2 



(A.18) 



5 The R amplitude has the phase of the 5-niatrix elements. We use the same conventions for the 
couplings as in m . 



14 



AJs, i, u) 



7T 2 VtU 2 



Q4g s d 



G 



Vt 2 



u- 



(A.19) 



where b\ and b 8 , in ( |A.18| , |A.19|) arise from the SM box diagram in Fig. lb. They are 
expressed in terms of the D functions defined in ED] and ET| as^ 



A = A(0, 0,0, m 2 H ,u, s; mf , mt, m*, m*) , (A. 20) 

Di = A(0, 0, 0, m 2 H , u, i; m t , rrit, m t , m t ) , ( A -21) 
Di = D i (0,0,0,m 2 H J,s;m t ,rnt,m t ,m t ) . (A. 22) 



h = m 2 (D + 2D 13 -2D 12 + D ) 

+t(-2D 133 + 2l! 122 - D 25 - 2D 23 - D 26 + 2% 2 + 2^ 24 - D l3 + 2^ 12 ) 
+u(-2D 123 + 2ll l23 - D 2A - 2D 26 - 2D 25 - D u + 2^ 26 

D 12 - 2D 13 - ]p Q - D 12 - i^ ) 

+s(-2D 233 + 2^ 223 - D 26 - 2D 23 - D 25 + 2^ 26 - D l3 + T) l2 ) 
-4D003 - 8D003 + 4^002 - 41)27 - AD 27 + AD 27 , (A.23) 



2 



8(D 2 33 - Diss - D 223 + D 23 - D 23 - D 26 ) + 2{D l3 - D 13 - D l2 ) . (A.24) 



and computed using the FF-package p8 



We have checked (numerically) that all gauge non-invariant terms arising from the indi- 
vidual diagrams in Fig. 3, cancel out when they are added. A further test is provided by 
the various Bose statistics relations (|A.15| - [A~T7D and ( |A.13j , |A.14|) , which are also perfectly 
satisfied. 



Appendix A3 : Kinematics for H+jet production in the parton 
model 

The basic parton model expression for the hadron-hadron collision AB — > Hf..., 
taking place through the subprocesses a + b — > H + / and illustrated in Fig. 2, is written 

cr(AB —> H f ...) = J2 dXadXb fa/A(x a )f b/B (x b )a(ab -> Hf) , (A.25) 

ab 

with / being a gluon or massless quark jet. Here f a /A(%a) is the distribution function of 
partons of type (a = g, q, q), in the hadron of type A. 

6 See in particular Eqs.(D.12) of @. 
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We next list a few formulae for the kinematics. The transverse momenta of the pro- 
duced H and the /-jet pt = Pth — Prf an d the energy E HT = \/p\ + m 2 H are described 
through 



._ . (3 T = p T /E TH = yjl-^jr, XfT = = PtXt ■ (A.26) 

The rapidities of the Higgs and the outgoing jet /, in the laboratory system are related to 
their energies and momenta along the beam-axis of hadron A, (taken as the i-axis) and 
to the corresponding production angles by 

e 2 - = l B+PHC066 a H , (A.27) 
E H -p H cos6 H 

e 2y f = E l± p l COs9 l 

E f - p f cos Of 

The center-of-mass rapidity y of the Hf pair, and the rapidities y* H , yj in cm., are defined 

as 

yH = y + y* H , Vf = y + y} • (A.29) 

The fractions of momenta carried by the incoming partons are given by 

*a = -£[e»" + PTe?t]=-=e* , (A.30) 
z 

= ^[e- yH +(3 T e- y f]=^=e-y , (A.31) 
Z ys 

while the Mandelstam invariants of the subprocesses satisfy 

s = M 2 = (p a + p b ) 2 = x a x b s = E 2 H [1 + (3 2 + 2p T cosh(Ay)] , (A.32) 

t = (ph - Pa) 2 = -E 2 TH f3 T {f3 T + e~ Ay ) , (A.33) 

u = {p f -Pa) 2 = -E 2 H [3 T ((3 T + e A y) , (A.34) 

§ 

t = - = x a x h , (A. 35) 
s 

where 

Ay = y H -y f = y* H - y) . (A.36) 
Denoting the cm. Higgs particle production angle by 9* and the cm. H velocity by 

P = Ph/K = , (A.37) 

T 1 1 I'll 
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we have also the expressions 



M = E TH (l+f3)coshy* H = p T (-+l)coshy} 



t 



u 



[m 



H 



-(l -cose*) 



1+COS0* 



Q * tanhy* H „ 
cos a = = — tanh y^ 



u — m 2 H 1 + (3 cos 9* 
i — m 2 H 1 — /? cos 0* 



The basic triple distribution is obtained from ( A. 25 ) to be 

da - rS% 
dp 2 T dy H dy f 13 



where, according to QA.25 ), is the total probability 

Sij = Sgg + Sgg + S ' gg + ^QQ 



with 



S'og = g{x a )g{xb) 



dcrjgg -» gg) 
di 



S 91 - ^ 



" vd^Ug ^ Hq) ,da(qg -> #g)' 

g[x a )q(x b ) ^ N(x„)5N' 



5, 



aY 



g{x a )q(x b ) da ^ 9q 2 H<1 ^ + q(x a )g(x b 
at 

q(x a )q(x b ) da qq 7 Hg + g(x a )g(x 6 ) 
at 



dt 




da(qg 




di 




da(qq — > 







where q = u,d, s, c, 6. From this basic distribution and imposing the cuts 

|yn| < Vff, |y/| < 17 > 



we get: 



The transverse energy distribution 

aV 



dxx 



1 I A M2XT C 

dyn I dy f — - — , 
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where M 2 is determined from ( |A.32| ) and the integration limits are 

P T x T 



Vfmin = max { In ( 2 _ x ^_ yn ) ; -Y f } , 

y fmax = min<!lnl — )'i Y f\, (A. 51) 



(3 T x T 



IjHrnax yHr 



(A,2, 



The rapidity distribution 

Since this is symmetric in yu, we may consider 



da 



J dxr J dyfM 2 xrSij , (A. 53) 



d\VH\ 

where the ^//--limits are given ( |A.51 ), while the limits of the xt integration are 

2tuh 



Lmin,exp 

o2 



XTmax = mm 

with 



5 " + i i 2e±VH ) (A. 54) 

I cosh y^' e2(i?±i/H)_i J v ' 



Ax = 4[e 2y ' + ^(1 - e 2 ( y /-^))] . (A.55) 
s 

The invariant mass distribution 

Using the cm. rapidity y defined above and ( |A.42|) , one obtains 

where the integration limits are 

r 1, Tr 1 (M 2 -xm 2 H \ , /V /s^ 
= ^(F g --ln X ; y / --l n ^ __ j; ln(-)|, 

iU. - max{-y,-iln X ; - y, (g=^) ; -]*($}, (A.57) 

r M 2 _ M^ + m 2 ^) . m|(l - e 2 ^/)) + ^ s 2yfl 
XmflX \m 2 H M±e~ 2Y f + sm 2 H ' 2M 2 ' M 2 

rm| M 4 e- 2y / + sm^ 2M 2 M 2 _ 2 y H ^ A , 

Xmm ~ ™*\M 2 'M 2 {s + m 2 H e- 2Y fym 2 H (l-e 2 ^+ Y f)) + ^ s 6 J (A j 
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where 



A 2 = m% ( e 2 ^ +y /) - l) 2 + 4M 4 e 2 ^ +y /) 



(A.59) 



The angular distribution is given by 



da 

dx 



(M 2 + m 2 H ) 



dM 2 I dy^—^Si,. 



(A.60) 



s(l + x) 2 v ' 

The y integration limits are as in (|A.57|), while for the M 2 integration we have the limits 



where 



M 



Mi 



mm 



in { X se 2Y »- S -e 2Y »- M 2 + ; M ' + 2 ;s} 



f Xm 2 H (e 2 ^ Y f +Y ^ - I) _ xm 2 H {e 2{Yf+YH) - I) 

I x 2 e 2(Y f +Y H ) _ i ' e 2(Y f +Y H ) _ ^2 



M 2 ; M' 2 



(A.61) 



M: 
A 3 



2Y,\2 



2/ 2 , Li, 

X {m H + se f 



Asm 2 H e 2Yf , 



M'l 



An 



2X 



2 2Yt\2 a 2 2 

m H + se ') — Asm H x e 



2 ^,2„2Y, 



f 



(A.62) 
(A.63) 
(A.64) 
(A.65) 
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Figure captions 



Fig.l Diagrams for top quark loop SM contributions to Hgg and Hggg couplings. 

Fig. 2 Diagrams for p + p — > H + jet + X. 

Fig. 3 Diagrams for gg — > Hg, gq — > f/g and — ► Hg. 

Fig.4 Rapidity distribution in p + p -> H + jet + X at LHC, for m H = 100 GeV . SM 
describes 1-loop SM predictions; SM e ff the large m t approximation to SM; +d G and — d G 
describe the Oqg contributions for d G = +10 -3 and d G = — 10 -3 respectively; and d G 
describes the O gg contribution for |g?g| = 10 -3 . 

Fig. 5 Angular distribution in p + p — > H + jet + X at LHC. See caption of Fig.4. 

Fig. 6 Transverse energy distribution in p+p — > H+jet + X at LHC. See caption of Fig.4. 

Fig. 7 Invariant mass distribution in p + p — > H + jet + X at LHC. See caption of Fig.4. 
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Fig. 4 



Rapidity distribution in p + p -> H + jet + X at LHC, for m H = 100 GeV. SM 
describes 1-loop SM predictions; SM e ff the large m t approximation to SM; +dc and 
-da describe the Oqg contributions for do = +10~ 3 and do = — 10~ 3 respectively; and 
dc describes the Oqg contribution for \da\ = 10~ 3 . 
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Fig. 6 

Transverse energy distribution in p + p — * H + jet + X at LHC. See caption of Fig.4. 
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